function [alpha L U] = solveAlpha (A,p)
%SOLVEALPHA Computes the alpha presented in C. M. Faria 2006.
%
%   INPUT:
%   A       :   Moment matrix.
%   p       :   Polinomial base at sample point.
%
%   OUTPUT:
%   alpha   :   Alpha solution
%   L,U     :   Lower and upper factorization of the .

    [L,U] = lu(A);
    alpha = (p/U) /L;
%    alpha = p/A;

end


